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A theory for light-induced current by strong optical pulses in molecular-tunneling junctions is de- 
scribed. We consider a molecular bridge represented by its highest occupied and lowest unoccupied 
levels, HOMO and LUMO, respectively. We take into account two types of couplings between the 
molecule and the metal leads: electron transfer that gives rise to net current in the biased junction 
and energy transfer between the molecule and electron-hole excitations in the leads. Using a Marko- 
vian approximation, we derive a closed system of equations for the expectation values of the relevant 
variables: populations and molecular polarization that are binary, and exciton populations that are 
tetradic in the annihilation and creation operators for electrons in the molecular states. We have 
proposed an optical control method using chirped pulses for enhancing charge transfer in unbiased 
junctions where the bridging molecule is characterized by a strong charge-transfer transition. An 
approximate analytical solution of the resulting dynamical equation is supported by a full numerical 
solution. When energy transfer between the molecule and electron-hole excitations in the leads is 
absent, the optical control problem for inducing charge transfer with linearly chirped pulse can be 
reduced to the Landau-Zener transition to a decaying level. When chirp is fast with respect to 
the rate of the electron transfer, the Landau theory is recovered. The proposed control mechanism 
is potentially useful for developing novel opto-electronic single-electron devices with optical gating 
based on molecular nanojunctions. 

PACS numbers: 73.63.Rt 



I. INTRODUCTION 

Molecular electronics is one of the most promising sub- 
stitutes for today's semiconductor electronics. In this 
relation molecular conduction nanojunctions have been 
under intense study in the last few years [l|, @, H, 0|- 
Recently, light induced switching behavior in the con- 
duction properties of molecular nanojunctions has been 
demonstrated @, 1, 0, 1, ©, 03 ■ 

However, the use of an external electromagnetic field 
as a controlling tool in the small nanogap between two 
metal leads is difficult to implement. Currently tech- 
niques available to achieve high spatial resolution with 
laser illumination are limited by diffraction to about half 
of the optical wavelength. The introduction of near-field 
scannin g op tical microscopes (NSOM) and tip-enhanced 
NSOM [ll| has extended the spatial resolution beyond 
the diffraction limit. The latter technique uses the 
strongly confined electromagnetic field generated by op- 
tically exciting surface plasmons localized at the apex 
of a sharp metallic tip, increasing spatial resolution to 
better than 10 nm [llj. Recently, spatial resolution at 
the atomic scale has been also achieved in the coupling 
of light to single molecules adsorbed on a surface, using 
scanning tunneling microscopy (STM) . 

If experimental setups that can couple biased molec- 
ular wires to the radiation field could be achieved, gen- 



eral questions concerning current through the molecu- 
lar nanojunctions in nonequilibrium situations come to 
mind. Recently Galperin and Nitzan investigated a class 
of molecules characterized by strong charge-transfer tran- 
sitions into their first excited state [l2j]. The dipole mo- 
ment of such molecules changes considerably upon exci- 
tation, expressing a strong shift of the electronic charge 
distribution. For example, the dipole moment of 4- 
dimethylamino-4'-nitrostilbene (DMEANS) is 7 D in the 
ground state and ~ 31 D in the first excited singlet state 
[13j. For all-trans retinal in polymethyl methacrylate 
films the dipole increases from ~6.6 to 19.8 D upon ex- 
citation to the 1 B U electronic state [l4| and 40 A CdSe 
nanocrystals change their dipole from ^0 to ~32 D upon 
excitation to their first excited state [l5| . In the indepen- 
dent electron picture this implies that either the highest 
occupied, or the lowest unoccupied, molecular orbitals 
(HOMO, |1), or LUMO, |2>, see FigHJ is dominated by 
atomic orbitals of larger amplitude (and better overlap 
with metal orbitals) on one side of the molecule than on 
the other and therefore, when used as molecular wires 
connecting two metal leads, stronger coupling to one of 
the leads. They have shown that when such molecular 
wire connects between two metal leads, weak steady-state 
optical pumping can create an internal driving force for 
charge flow between the leads. 

A theory of light-induced effects by weak CW radiation 
in molecular conduction was developed in Ref . [la ] . How- 
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ever, there are reasons to consider also molecular junc- 
tions subjected to strong electromagnetic fields. First, 
the structure of such junctions is compatible with config- 
urations considered for large electromagnetic field as in 
tip enhanced NSOM [ll| . Secondly, it was demonstrated 
in Ref. [Il[ that the combination of near-field optics and 
ultrafast spectroscopy is readily achieved, and the obser- 
vation of photo-induced processes, such as charge trans- 
fer, energy transfer or isomerization reactions on the 
nanoscale is feasible [l7j- Third, consideration of junc- 
tion stability and integrity suggests that strong radiation 
fields should be applied as sequences of well separated 
pulses to allow for sufficient relaxation and heat dissi- 
pation. Finally, consideration of strong time dependent 
pulses makes it possible to study ways to optimize the de- 
sired effect, here the light induced electron tunneling, i.e. 
to explore possibilities of coherent control of charge flow 
between the leads. Our objective in the present work 
is to extend the theory of Refs.fH, [l6[ to strong fields 
and to apply the theory to studies of coherent control of 
nanojunction transport. 

While these problems are of general and fundamental 
interest, we note that this study is related to efforts to 
develop novel optoelectronic single-electron devices, such 
as a photon-electron conversion device, optical memory, 
and single-electron transistors with optical gating [|| . In 
addition, the potential significance of molecular nano- 
junctions for device applications lies in the possibility 
of creating all-optical switches [I8| that could be incor- 
porated in future generations of optical communications 
systems. It is conceivable that these devices will employ 
coherent optical manipulations, because the speed of co- 
herent manipulations greatly exceeds that of currently 
available electronic devices. 

The outline of the paper is as follows. In Sec|TT] we 
introduce our model. In Sec lIIII we derive a closed set 
of equations for the expectation values of binary and 
tetradic variables of the annihilation and creation opera- 
tors for electrons in molecular states |1) and |2), and get 
formulas for the current and charge transferred during 
the electromagnetic pulse action. In Sec lIVI we calculate 
a current induced by quasistationary intense light pulse. 
Optical control of current and transferred charge with 
chirped pulses is considered in Sec[V] We summarize our 
results in Sec lVIl In the Appendices we show that in the 
absence of the radiative and nonradiative energy transfer 
couplings, the equations of motion derived in the paper 
lead to the well known Landauer formula for the current 
and present auxiliary calculations. 



II. THE MODEL HAMILTONIAN 

We consider a system that comprises a molecule 
represented by its highest occupied molecular orbital 
(HOMO), |1) , and lowest unoccupied molecular orbital 
(LUMO), 1 2) , positioned between two leads represented 
by free electron reservoirs L and R and interacting with 




FIG. 1: A model for light induced effects in molecular con- 
duction. The right (R = \{r}) ) and left (L = |{/}>) manifolds 
represent the two metal leads characterized by electrochem- 
ical potentials ^ir and \il respectively. The molecule is rep- 
resented by its highest occupied molecular orbital (HOMO), 
|1), and lowest unoccupied molecular orbital (LUMO), |2). 



the radiation field (FigQJ. In the independent elec- 
tron picture a transition between the ground and excited 
molecular states corresponds to transfer of an electron 
between levels |1) and |2). The electron reservoirs (leads) 
are characterized by their electronic chemical potentials 
/il and where the difference /i£ — [Xr = e$ is the 
imposed voltage bias. The Hamiltonian is 

H = H + V (1) 

where 

m=l,2 ke{L,R} 

contains additivcly terms that correspond to the isolated 
molecule (m) and the free leads (k). Here hi — cf£i is 
the population operator in state i, the operators c and 
c + are annihilation and creation operators of an electron 
in the various states. 

The interaction term V can be written as 

V = V M + Vp + Vn (3) 

where Vp accounts for the effect of the external radiation 
field. The latter is represented by the (classical) function 

E(r,t) = E (+ )(i)+E (-) (i) = -e£(t) exp[~iujt+i(p(t)}+c.c. 

(4) 

characterized by the pulse envelope £(t), carrier fre- 
quency u> and (possibly) time dependent phase <p{t). The 
time dependent phase corresponds to time evolution of 
the pulse frequency (chirp) u>(t) = to — dip(t)/dt. Intro- 
ducing bilinear operators of the excitonic type 
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Hi = £ tci, b M = c+c 2 , b+ = c+Cj = b jh (i ^ j), (5) 

the molecule-radiation field coupling, Vp, can be written 
as follows in the resonance or rotating wave approxima- 
tion (RWA) 

V P = --(d ■ e){6+ £{t) exp[-iu>t + i<p(t)] + h.c.} (6) 

where d is the transition dipole moment. 

The other terms in Eq.([3]) describe coupling between 
the molecule and the metal electronic subsystems. In 
terms of the excitonic operators defined in Eq.([5]), they 
are given by 

Vm= £ £ (V^ K h mk + h.c), (7) 

K=L,B.m=l,2;keK 

Vn= E E (V^h^+V^bubi,,), (8) 

L and R denote the left and right leads, respectively, and 
h.c. denotes Hermitian conjugate. Vm and Vn, Eqs. (|7|) 
and (O, respectively, denote two types of couplings be- 
tween the molecule and the metal leads: Vm describes 
electron transfer that gives rise to net current in the 
biased junction, while Vn describes energy transfer be- 
tween the molecule and electron-hole excitations in the 
leads. The latter interaction strongly affects the lifetime 
of excited molecules near metal surfaces. Vn is written 
in the near field approximation, disregarding retardation 
effects that will be important at large molecule-lead dis- 
tances. 



III. EQUATIONS OF MOTION 

The physics of the system can be described within 
different approaches. One is the method of nonequi- 



librium Green's functions [H, Ell IHl- It has advan- 
tages of a formal treatment due to the possibility of a 
diagrammatic representation, and it is particularly well 
suited for stationary processes where the Dyson equa- 
tion can be cast in the energy representation. For time- 
dependent processes, such as are the subject of this work, 
a method based on the equations of motion for the expec- 
tation values of the operators provides a more transpar- 
ent approach, since the quantities are more directly re- 
lated to physical observables. Such a method is adopted 
here. Using a Markovian approximation for the relax- 
ation induced by the molecule-metal leads coupling, we 
derive a closed set of equations for the expectation val- 
ues of binary (n m ) = n m and (6a/) = pm, and tetradic 
(b M bM) = Nm variables of the annihilation and creation 
operators for electrons in molecular states |I) and |2). 
The first expression is simply the population of electrons 
in molecular state m, the second gives the molecular po- 
larization and the third represents the molecular excita- 
tion, referred to below as the molecular exciton popula- 
tion. 

Using the Heisenberg equations of motion one obtains 
the equation for the expectation value of any operator F 

j t (F) = ~{[H + V,F]) = ~Tr{[H Q + V,F]p) (9) 

where p is the density matrix. Straightforward operator 
algebra manipulations yield for n m and p m in RWA 



= {-l) m Im{n* (t)p M exp[iujt-iLp{t)}}- 

E E^T } vW- 

K=L,Rk£K 

-l Im E E [^mV^ K \b M bi, k ) + 
K=L,Rk^k'eK 

+ SlrnV^ K) (b M b+ k y] (10) 



dpM 
dt 



\[e x - e 2 )p M + |Q(t)exp[-icjt + i^(t)](rii - n 2 ) + j- J2k=l,rJ2 

+ i T,K=L,p.T, k ^k'eK V kk> ] (bk'k(n 2 - fix)} 



( y(MK) (b+) 
k£K\ v kl \ u k2l 



V, 



(MK) 



2k 



(hi)) 



(11) 



where ft(t) = (d • e)£(t)/H is the Rabi frequency. The 
equations of motion include couplings to additional cor- 
relations of the second order (b m k) due to the electron- 
transfer interaction Vm, and to higher-order correlations 
(&M&fc' fc ) etc. due to the energy transfer Vn- To obtain 



expressions for these correlations, we now compute their 
equations of motion, using the Markovian approximation 
for the relaxations induced by the molecule-metal leads 
couplings, Vm and Vn- In this work we assume that the 
relaxation processes due to Vm and Vn are not inter- 
dependent and also do not depend on the external elec- 
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tromagnetic field. We shall discuss the last approach in 
SecEI 



becomes 

( bmk ) = {bkm y = U J2 v^l K) (ci,c m )(t)x 



A. Calculation of terms associated with the 
electron transfer interaction Vm in the equations for 

and pm 

In evaluating the effect of the relaxation processes as- 
sociated with the electron transfer and energy trans- 
fer terms in the Hamiltonian, Vm and Vjv respectively, 
we make the approximation (known as the non-crossing 
approximation) that these processes do not affect each 
other. A similar assumption is made with respect to the 
effect of the external field. With this in mind we consider 
the expectation values (bkm) and (b m k) that enter the 
terms containing Vm on the right-hand side of Eqs. (fT0|) 
and and omit Vp and Vm terms in the equations of 
motion that describe their evolution. This leads to 



Jt {bmk) = li 



e k -e m ){b mk ) + - V$ 



k \ c m'°r, 



l v (MK) f , v 



(12) 



where we assumed that the leads are in equilibrium with 
the expectation values 



(13) 



ihP 



+ tvK5(e k - e m >)}- 



VT' W)[ 



ihP 



£k ~ £ r . 



+ hir6(e k - e m )}} (15) 



where P denotes the principal value. Substituting the 
last result into the corresponding terms containing Vm 
on the right-hand side of Eqs. p!)]) and (fTTj) and keeping 
only resonant terms, we have 



Im £ J2 V krn(b mk )^ E [n m (t)T M ,, 

K=L,R 



K=L,RkeK 



W M K,; 



(16) 



and 



K=L,Rk£K 



-PM 



(*) E \\{ T MK,l+T MK i)+i&MK] (17) 



K=L,R 



where 



MK,i 



2?r 



kGK 



(MK), 2 



S(e k 



(18) 



f K (e k ) = [exp((e fc - fi K )/k B T) + l]" 1 

is the Fermi function, S k k' is the Kronecker delta. For- 
mally integrating Eq. (fT2"]) . we get 



(bmk) = T dt' exp[-(£fc - e m )(t - t')]x 

x[ e ^r ) (c+,c ro) (o-^r ) /K( £fe )] 



(14) 



In the absence of the Vm coupling this results in a 
set of integro-differential equations for (c+c m ) = n m 
and (£^£2) = pm- The dynamics contains mem- 
ory effects, and therefore is non-Markovian. Next 
we make a Markovian approximation by transforming 



to the interaction representation: (c^/C m )(t') 



*MK 



IV, 



(Mi?) ,2 



fel 



w. 



(MK) |2 



21, 



k£K 



£k - ei 



£k ~ £2 



(19) 



is the correction to the frequency of molecular transition 
(e% — £\)/h due to electron transfer between the molecule 
and lead K, 



W M k,, 



2tt 



sEW 6 *)!^! 2 ^-^) (2°) 



keK 
fK\£m)^MK 



(21) 



'it') expire, 



It should be noted that the second equality, Eq. (|21jl . is 
valid only provided that molecular state e m is far from 
the Fermi level of lead K. The point is that the position 
of e m influences on the bath correlation frequency u> c , 
for the relaxation parameter WMK,m (^c is the frequency 
interval at which the interaction of a system (molecule) 
with bath (metal leads) is essentially changed) . Really, if 
molecular states e m are far from the Fermi levels of both 
leads, u> c for WMK,m is the same as that for TMK,m, 
which is determined by the frequency interval for the 



e m )t'] and assuming that system-bath interaction matrix element V k ^ K ^ and the 



slowly varying function (c^ l ,c m ) mt (t / ) can be moved as 
(c m ,c m ) int (t) to outside the integral 0]. Ea.([TlJ) then 



density of states of metal leads. For this case uj c can be 
evaluated as 1-10 eV 1201 . The situation is different if we 



5 



assume that the molecular level position is pinned to the 
Fermi energy of a lead. In the latter case uj c for Wmk,™, 
is determined also by the frequency interval at which 
is essentially changed that is <~ fc^T/ft = 0.026 
eV for room temperature (see Eq. (f20"]) ). Since the value 
of lu c places a limit on the used approximation, accord- 
ing to which the relaxation parameters do not depend 
on exciting electromagnetic field (see Sec lVI[) . one can 
use Eq. pTj) only in the case when the bath correlation 
frequency uj c is the same for both WuK,m and ^MK,m ■ 
One can easily sec from Eqs.lfTD]). (fTTjl . (Ti!)|) and (fT?) 
that in the absence of energy transfer (V/v)> equations for 
the populations of molecular states and molecular polar- 
ization form a closed set of the equations of motion. 

B. Calculation of terms related to energy transfer 
in the equations for n m and pm 

The calculation of terms related to energy transfer in 
Eqs.|T0l) and ([TT]) is similar to that of Sec lIII Al Invok- 
ing again the non-crossing approximation by omitting Vp 
and Vm terms in the equations of motion for the expec- 
tation values (&M^fe/fc) and (b k ' k {n 2 — ni)), which appear 
on the right-hand sides of Eqs. (fT0|) and (flTj) . respectively, 
we get 

d i 

+ ^ K) {f K (e k )[l .f K (sk')}{b + M b M )- 
- f K (e k ')[l- fK{e k )]{b M b + M )} (22) 



-^(bk'k{n 2 - fix)) = -(e k - e k ')(b k ' k (h 2 - hi)) + 

+ ^ N k K) {f K (e k/ )[l-f K (s k )} + 
+ f K {e k )[l- f K {e k -)]}VM (23) 



Formally integrating the last equations, performing 
Markovian approximation and substituting the results 
into the corresponding terms containing Vn on the right- 
hand side of Eqs. fTO]) and (|lip. we obtain 



E E v.T'ibMbU) 

K=L,Rk^k'eK 

= ^2 [B NK {ei- e 2 ,^K){bMb\ I ) 

K=L,R 

-B NK (s 2 -e 1 , f i K )(b M b M )] (24) 



E 



E 1 

K=L,Rk^k'eK 



V kk' K ( b k'k{n 2 - fix)) = -p M (t) E i iA NK + -[B NK {ei - £2,(i>k) + B NK (e 2 - El,flK)]} (25) 



K=L,R 



where 



2tt 



E w£ 



(JVif) |2 , 



k^k'eK 

x5{e k - e k - + e m - e n )f K (e k )[l - /#(£#)] 



(26) 



and 



^nk = \p E {fK(e k ')[l- fi<(e k )} + 



k^k'£K 

+ f K (e k )[l- f K (e k <)}} 



£k - £k> + £ 2 - El 



(27) 



is the correction to the frequency of molecular transition 
(e 2 — e±)/h due to energy transfer between the molecule 



and lead K. In deriving Eqs. (l24]) and ([25]) we have used 
the arguments, which are similar to those used above in 
the derivation of Eqs.flTBJ) and (fT7| . 



One can see from Eqs. (fTU|) and (J2U) that in the presence 
of energy transfer (Vn), equations for the populations of 
molecular states and molecular polarization do not form a 
closed set of the equations of motion. They must be sup- 
plemented, at least, with equations for the expectation 
values of tetradic variables (&m&m) an d (^m^m) — Nm, 
where (b\ib M ) and Nm are related by the following equa- 
tion: 



(b M b M ) = N M -n 2 + n 1 



(28) 
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C. Equation for (&Af&jvf) 

Using Eq. ([9]) , straightforward operator algebra manip- 
ulations yield for (&m^m) m RWA 



d(b M b 



Ml 



dt 



K=L,R k£K 

-vL MK \b klb + M ))+ 

+ E E v£ K \b M bt, k )- 

K=L,Rk^k'eK 

— Im{f2*(i) exp[iwf — i(p(t)]pM} 



(29) 



where the second term on the right-hand side of Eq. (j29|) 
has been calculated above, Eq. (|24j) . The first term on the 
right-hand sides of Eq. (|2"9"|) is associated with the electron 
transfer process. To evaluate it we consider the equa- 
tions of motion for the expectation values (&2fe&i) and 
(bkib M ), omitting Vp and Vm interactions and keeping 
only resonant terms: 



d(b 2k b 



•2kO M ) _ 



= Ue k - e 1 )(b 2k b+) + ^ K) [(b M b+)- 



dt h 

-f K (Sk){l 



- n 2 ) 



(30) 



d(bkib 



ill 



w/ =\{s2-e k ){b kl bl II 



-IV^ K) [(1- f K (e k )) ni - (b M b M )] 



(31) 



Integrating Eqs. (|30f and (f3T|) . performing Markovian ap- 
proximation and substituting the results into the first 
term on the right-hand side of Eq. ([29| , we get 



-| Im E J2(vL MK) (h k b + M )-vi MK \b kl bt I )) = 

K=L,Rk£K 

= E i-( b M b tl)( r MKA + T M K,2) + (1 - n 2 )W M K,l + 
K=L,R 

+ {Tmk,2 - W MK ,2)ni] (32) 



where T M K,m and Wmk. 

m were defined in Eqs. (|18p . (120|) 

and (ED. 



D. A closed set of the equations of motion 



using Eqs.dujj, cui, dinD, cud, ©, cu), qd, 

(HI, (123, (EH]), (HZ]), (EHD, ©, (EH and switching to the 
system that rotates with instantaneous frequency w(t), 
PAi(t) = pm(t) exp{i[ujt — f(t)]}, we obtain a closed set 
of equations for the quantities that vary slowly with time 
during the period of a light wave 



dn Tt 



dt 



Ml f W Mm - TumUm - {-l) m [B N {e 2 - ex)N M - {n x - n 2 + N M )B N (e 1 - e 2 )] (33) 

I 



=i[u(t)-u ]pM+~ft(t){ni-n 2 )-^r MN p M (34) 



dN 



dt 



= (Tmi - W M x)n 2 + Wmh(1 - ni) + Im[ft*(t)p M ]- 



— [Bat(£2 — £l) + Tmi + rAf 2 ]A^M- 

+ ("i - n 2 + N M )B N {ei - e 2 ) 



where 



K=L,R 



K=L,R 



K=L,R 



(35) 



TMm — ^ rjvfK.mj I^Mm — ^ WuK,m (36) 



^ BNK( £ m — E ni^K) — B N (e m — £ n ) (37) 



Tmw — Tmi + Tm2 + Bn(e 2 — £i) + Bjsr(si — e 2 ) (38) 



and w = (^2 - £i)/ft + Y.k=l,r^ nk + ^mk) is the 
frequency of the molecular transition with the corrections 
due to energy and electron transfer between the molecule 
and the leads. 

As indicated above, equations for the populations of 
molecular states and molecular polarization, Eqs. (|33p 
and (|34|) . form a closed set of the equations of motion 
if the energy transfer is absent (B]y(e m — e n ) = 0). 
When energy transfer is present they must be supple- 
mented with Eq. (|35|) for the exciton population. On the 
other hand, in the absence of electron transfer (WMm — 
^Mm = 0) Eq. ([33|) coincides with Eq. (f3"3"|) for n 2 , which 
implies that Nm = n 2 . Indeed, Nm = (c 2 £icfc 2 ) — 
(h 2 (l — n\)) = (h 2 ) = (h 2 ) = n 2 when the electron popu- 
lation on the molecule is conserved, i.e. when Vm = 0. It 
is the combined effect of the electron and energy transfer, 
represented by the terms Vm and Vm in the Hamiltonian, 
that leads to the need to include Eq. (|3"B"|) in the closed set 
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of the equations of motion. 



^ = Im[fi*(i)p M ] - T M2 n 2 - B N (e 2 - £l )JV M (42) 



E. Calculation of current and transferred charge 

The electronic current I is given by the rate at which 
the number of electrons changes in any of leads, e.g. fl9l . 




(39) 



Evaluating the commutator in Eq.([39|). we get 



2e 

1 



E E^vW 

m=l,2 fcei 



E 



' > lnm(i)rMl,m - WmL, 
m=l,2 



(40) 



where we used Eq. dlGp . Correspondingly, the charge 
transferred during an electromagnetic pulse of finite du- 
ration is given by Q = I{t)dt. 

In the Appendix A we show that in the absence of the 
radiative and nonradiative energy transfer couplings, Vp 
and VW, Eqs.dHJ), (HJ), ©, ([36]) and gO]) lead to 
the well known Landauer formula for the current [2l| . 



IV. CURRENT INDUCED BY A 
QUASISTATIONARY LIGHT PULSE 

In this section we calculate the current induced in 
a molecular nanojunctions by a strong quasistationary 
light pulse. Here and in the next section we assume 
that the molecular energy gap e 2 — £i is much larger 
than the voltage bias hl — Hr and that the HOMO 
and LUMO energies, e± and e 2 , are positioned rather 
far kpT) from the Fermi levels of both leads, so 
that the dark (Landauer) current through the junction 
is small and may be disregarded. Using for this situation 
Wmk,i ~ T M k,i, Wmk.2 ~ B NK (ei - s 2 ,fi K ) ~ 0, we 
obtain from Eqs.([33|), (jM )) . (|33|) and P0f 



drii 
~dT 



Im[n*(t)p M ] + r M i(l - n x ) + B N (e 2 



(41) 



= i[u)(t)-u)o]pM+^{t){nx-n 2 )-^MNPM (43) 



dN M 
dt 



lm[n*(t)p M ] ~ [B N (e 2 - ei) + T M i + T M 2]N M 

(44) 



/ = e[(m - l)T ML> i + n 2 T ML ^ 2 ] 



(45) 



One can see from Eq. (l4"5]) that the current strongly in- 
creases when n 2 , 1 — ri\ ~ 1, which can be realized for 
strong light fields. If we further assume that the pulse 
amplitude £{t) and frequency co(t) change slowly on the 
time scale of all relaxation times as well as the reciprocal 
Rabi frequency, one can put all time derivatives on the 
left-hand sides of Eqs.gT]), (gU), (03]) and (gU) equal to 
zero, and the resulting stationary equations can be easily 
solved 



"2 



o 2 (t)(r M2 



Ml 



)/(4T 



.1/2 J 



m = 1 - n% 



M2 



Ml 



N. 



A/2 



M 



-n 2 



Ml 



M2 



(46) 



(47) 



(48) 



pm 



fi(f) 



i(r M jv/ 2 ) - [u(t) - w ] 



n 2 (t) 



(49) 

This solution corresponds to the molecular level and exci- 
ton populations as well as polarization adiabatically fol- 
lowing the optical pulse. Substituting Eas.(|3S|) and (j4"Tl) 
into Eq. (|4"S"]) . we get 



I{t) 



Tm2 + TMI ^ 2 (t)(^ML,2^MR.l — ^ MLS^ MR.2.) 



4T M ir M 2 n 2 (t) 



4r A/1 rj. 



+ (^) 2 + [v - cu(t)Y 



(50) 



At steady-state (ui(t) = u>, il(t) — fl) and small fields, 
2 iry+r«o 2 (T MN /2) 2 , this becomes 

41 Mil M2 v ' ' ' 
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/ = 



efi 2 



M2 



Ml 



^ML,2^MR,1 



ML,1^MR,2 



4 (T MN /2) 2 + (uj - ujf 



TmxT 



A12 



(51) 



The last equation is similar to Eq.(50) of Ref.[lj| with 
the only difference that the latter corresponds to the 
substitution of the sum Tm2 + Ta/i on the right-hand 
side of Eq.([5U by T MN = T M2 + T m + B N (e 2 - e ± ) > 
Tm2 + rji/i- The difference may arise from the fact that 
Eq.(50) of Ref.[l6[ is obtained in the much used approxi- 
mation of strong dephasing [22] Nm = (^2(1 — fti)) ~ 
"2(1 — Ti%). For small fields the latter term is of or- 
der fl 4 , since n 2 , 1 — n\ ~ il 2 . As a matter of fact, 
when the exciting field is weak, one can neglect the term 
Bn(£2~ £1)^2(1-^1) ~ ^ 4 with respect to Fmi(1— n\) ~ 
O 2 and Tm2T12 ~ ^ 2 on the right-hand sides of Eqs. (|4"Tj) 



and (|42p . respectively. In other words, the approxima- 
tion of strong dephasing Nm w ^2(1 ~ n i) disregards the 
depletion of state 2 due to energy transfer for small fields, 
and therefore results in some overestimation of the cur- 
rent. In contrast, our present approach takes the tetradic 
variable Nm into account exactly (in the framework of 
the Markovian approximation) and does describe the de- 
pletion of state 2 due to energy transfer. 

For strong fields and near resonance excitation, 

^(t) ( Sf » (r MW /2) 2 ,N-^(t)] 2 , Eqs.®, 
(|47|) . (14511 and ([49]) describe the saturation effect 



r 



m=n 2 = 



Ml 



netic field limit. As mentioned in the introduction, future 
generations of optical communication systems will em- 
ploy coherent optical manipulations whose speed greatly 
exceeds that of currently available electronic devices. We 
next consider such coherent control processes. 

Two well-known procedures based on a coherent ex- 
citation can, in principle, produce complete population 
inversion in an ensemble of two-level atoms. One of them 
is the 7r-pulse excitation [231 ] . which makes use of the 
Rabi population oscillations. This approach has been 
successfully demonstrated in atoms as well as semicon- 
ductor quan tum dots, often referred to as artificial atoms 
[2H HHliE H3- The main disadvantage of the 7r-pulse 
excitation method is the requirement for resonant light 
source and the need of precise control of the pulse area 

The second procedure, known as adiabatic rapid pas- 
sage ( ARP) [13, [H, [H, M, HE H IH ll , enables us to 
transfer the entire population from ground |1 > to the 
excited \2 > electronic state. It is based on sweeping the 
pulse frequency through a resonance. The mechanism 
of ARP can be explained by avoided crossing of dressed 
(adiabatic) states 



M2 



Ml 



(52) 



N 



FmiT 



M2 



M 



(fM2 + Tmi) 5 



Pm 



2TmiTm2 i(T M n/2) - [w(t) - lo ] 



(Tm2 + fjlfl) 2 



n(t) 



and Eq. ([40[) gives 



7 = e 



^ML,2^MR,1 — FmL,iFmR,2 



M2 



+ r 



Ml 



(53) 



(54) 



(55) 



Eas. ([5D|) and Eg. ([55)1 show that the optically induced cur- 
rent increases linearly with the pulse intensity ~ f2 2 for 
weak fields, and saturates at the maximal value given by 
Eq. ([55|) . for strong fields. As is easy to see from Eqs. ([52"|) 
and ([53")) . in the latter case Nm = n 2 (l — n\), since the 
strong dephasing limit where Pm can be disregarded is 
realized under saturation effect (see Eq.{54])). 



V. OPTICAL CONTROL OF CURRENT AND 
TRANSFERRED CHARGE WITH CHIRPED 
PULSES 



In the previous section we have generalized the results 
of Refs.[l2,[l6| to the quasistationary strong electromag- 



net) = Bini?(i)|l > + costf(t)|2 > 
$_(*) = costf(*)|l > -sin<?(t)|2 > 



(56) 



as a function of the instantaneous laser pulse frequency 
ui(t) [23]. Here the mixing angle $(t) is defined (modulo 
71") as d — (1/2) arctan ^ n °^ t ) , where Q(t) is the Rabi 
frequency. During the excitation, the mixing angle ro- 
tates clockwise from $(— 00) = tt/2 to i3(+oo) = and 
the composition of adiabatic states changes accordingly. 
In particular, starting from state |1 >, the system follows 
the adiabatic (dressed) state ■£+(£) and eventually ends 
up in state |2 > [3l|. A scheme based on ARP is robust 
since it is insensitive to pulse area and the precise loca- 
tion of the resonance. Therefore, we shall focus in what 
follows on ARP as a way to control optically induced 
charge transfer in molecular nano junctions. The appli- 
cation of our formalism to the coherent optoelectronic 
properties of nano junctions with quantum dots, using tt- 
pulse excitation, [23], will be analyzed elsewhere. 

As a particular example we shall consider a light- 
induced charge transfer in molecular nanojunctions, us- 
ing linear chirped pulses uj(t) — u> — p,(t — to) where 
p, = d 2 ip(t)/dt 2 = const. 
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FIG. 2: The charge Q transferred after the completion of 
the pulse action as a function of the linear chirp rate fi' = 
fi/uio (u>o is defined below Eq. (|38|l 'l. Other parameters are 
as follows u = ujq = 3 eV, Vm%/^o = 0.03, Tmi/uq = 0.04, 

r MN = r M1 +r M2 , d£o/hu> = 0.2, 5/^0 = 0.1, r M1 , : = 0.01 

eV, r ML , 2 = 0.02 eV. 



A. Numerical results 

Figures and [3] show the influence of fi, the chirp rate 
in the time domain, on the charge transferred during one 
electromagnetic pulse action. These results are obtained 
by numerical solution of Eqs.flH]), ([32]). (ITS)) . (|3I|1 and 
(|4"S")) for a Gaussian pulse of the shape 

E{t) = E (t) exp (»> (*)) = £ exp[-i(<5 2 - ip)(t - t ) 2 ] 

(57) 

and are displayed as a function of \i. We see that pulse 
chirping can increase the transferred charge (Fig|5]) and 
the induced current (Fig[3]) that can be explained by sig- 
natures of ARP (see below). 

If chirped pulses are obtained by changing the separa- 
tion of pulse compression gratings, the parameters 5 and 
jl are determined by the formulae [35l l.36| : 

S 2 = 2r 2 [r^ + W' 2 (u)}-\ 
/2 = -4$"(o;)[r p 4 + 4$" 2 (a;)]- 1 ) (58) 
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FIG. 3: The current 7 as a function of time r = ujot for the 
linear chirp rate Ji/ul = (A), 0.07 (B) and 0.15 (C). Other 
parameters are identical to those of FigGD The figure illus- 
trates how signatures of ARP increase the induced current. 
Inset, the square of electric field amplitude of the exciting 
pulse in arbitrary units. 



where r p o = t p o / V21n2, t p o is the pulse duration of 
the corresponding transform-limited pulse, and $" (lo) 
is the chirp rate in the frequency domain. The latter 
is defined by writing the electric field at frequency ui 
as exp [i$ ((D)], and expanding phase term <E»(a>) in 

a Taylor series about the carrier frequency ui = 
$(w) + (l/2)$"(w)(<D - uj) 2 + ... Note that the local field 
in the junction reflects also plasmon excitation in the 
leads, and taking the incident pulse shape as affected 
only by the compression gratings used disregards the pos- 
sible contribution of the near field responce of plasmonic 
excitations in the leads [ij], H3 • Such effects will be con- 
sidered elsewhere. 

Figs|l]and[5]show the calculation results of the trans- 
ferred charge Q as a function of the chirp rate in the fre- 
quency domain $" (v) = 4ir 2 &" (ui). The calculated de- 
pendences Q(Q" (v)) for curves A, B and C are confined 
to the values of an argument |$" (v)\ > corresponding 
to dSo/fkuQ < 0.3 (d is the molecular dipole moment, cf 
Eq.©), since our theory uses RWA. One can see that Q 
grows rapidly for small |$" (y)\ . The growth of Q slows 
down for moderate |4>" (z/)| , and then Q tends to a con- 
stant value for large | $" {v) \ . The larger is pulse energy, 
the larger is the value of |$" {v)\, at which the growth of 
Q slows down. FigQ] corresponds to the absence of the 
energy transfer (Bjv(£2 — £i) = 0), and Fig. [S] illustrates 
the influence of the energy transfer (Bn{s2 — £%) 7^ 0), 
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FIG. 4: The charge transferred after the completion of the 
pulse action as a function of the chirp rate in the frequency 
domain $" (v). The value of d£o/hu = 0.6 (A), 0.5 (B), 
0.4 (C) and 0.3 (D) for the transform-limited pulse. In 
the course of chirping the pulse energy is conserved so that 



2 L'pO ' 



creases when |$" (y) \ increases; r p o 
are identical to those of Fig[2] 



= const, and £o de- 
ll fs. Other parameters 



which diminishes the corresponding values of Q (sec also 
Figl5]). The behavior and values of Q shown in FigsHHU 
[5] and [6] can be rationalized by the theoretical considera- 
tion below. 

Fig. [7] illustrates the influence of detuning between the 
carrier pulse frequency ui and the corrected frequency of 
the molecular transition loq on the transferred charge Q. 

To end this section we note that the current that cor- 
responds to the expectation value of Q = 0.5 • 10~ 19 C 
per pulse (corresponding to curve D in Fig|4]) and to an 
estimated pulse repetition frequency of 82 MHz [27| re- 
sults in a small but measurable value of about 4 x 10~ 12 
ampere. 



FIG. 5: The charge transferred after the completion of the 
pulse action as a function of the chirp rate in the frequency 
domain <E>" (v) in the presence of energy transfer Bn(£2 — 
£i)/wo = 0.01. Other parameters are identical to those of 
Figgl 
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B. Analytical consideration 

The problem under consideration above can be solved 
analytically in certain conditions. 



FIG. 6: The charge transferred after the completion of the 
pulse action as a function of the chirp rate in the frequency 
domain $" (v) when d£o/hwo = 0.1 in the absence (Bn(e2 — 
£i)/wo = 0, curve A) and presence (Bjv(e2 — Ei)/uo = 0.01, 
curve B) of energy transfer. Other parameters are identical 
to those of Fig|4] 



1. Chirped pulse control of charge transfer in molecular 
nanojunctions as the Landau-Zener transition to a decaying 

level 

Consider first an excitation of the molecular nanojunc- 
tion with a linear chirped pulse uj(t) = u>o — jit (to = 0, 
lu = loq) of a constant amplitude = const) in the 
absence of energy transfer: Bjy(e2 — £i) = Ank = 0, 



Tmjv = T M 2 + Tmi (rjvfiv was defined by Eq.JSS])). If in 
addition, Tmi = Ta/2 = ^m an d provided that level 1 is 
below and level 2 is above both Fermi energies, then it 
can be shown that m = 1 — ni (see Appendix B) and 

/oo 
n 2 (t)dt (59) 
OO 
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FIG. 7: Influence of the frequency detuning uoq — uo on the 
charge transferred after the completion of the pulse action for 
d£o/hcuo = 0.1. (u>o — ^)/<^o = (solid line) and 0.05 (dashed 
line). Other parameters are identical to those of Fig|4] 



Under these conditions our electron problem (Eqs. (j¥T| )- 
(|44l) ) becomes mathematically equivalent to the Landau- 
Zener transition to a decaying level[45| solved analyt- 
ically by Akulin and Schlcich [38]. The magnitude 
-Tml,2 n 2 (t)dt on the right-hand side of Eg. ([55)1 rep- 
resents the expectation value of the number of electrons 
passed from the molecule to the left lead after the com- 
pletion of the pulse action, and Tml,i J_ oci n 2 (t)dt = 

^ML,i (1 — ni(t))dt is the same for the electrons 
passed from the left lead to the molecule. 

Using Eq.® and Eq.(25) of Ref.[H for the magni- 
tude I as = Tju J™ n,2(t)dt, we obtain in terms of our 
representation 



(air)^ 1 sinh(7ra). The expression in the brackets on the 
right-hand side of Eq. (j6"Tj) is simply the probability of 
Landau-Zener transition, which indeed is identical to 
the probability of the electron transfer from the excited 
molecule to the leads in the case of fast passage through 
the resonance. Indeed, in this case Tm J_ ao n 2(t)dt = 

Tm / °° ^2(0) exp(— Titftydt = n 2 (0) where n 2 (0) is the 
population of molecular state 2 immediately following the 
passage through the resonance. The highest charge trans- 
fer is therefore obtained if 712(0) = 1. Eq. (f6~lj) shows that 
712(0) approaches 1 for strong interaction, ttVL 2 ^> 2|/i|, 
which corresponds to adiabatic rapid passage (ARP). In 
other words, when the interaction with light is short in 
comparison with the electron transfer, the transferred 
charge is maximal when ARP is realized. Really, Q — > |e| 

Jj? rA/z, 7 2— Fm L.l 



1. This issue is of importance for de- 
veloping single-electron devices with optical gating based 
on molecular nanojunctions. 



When 



A 

Tp7 



1, the magnitude I as is given by [36 



Ias 



7rrj\/f2 



2|jS| ^m + W 



and we get a simple formula for the charge transferred in 
the course of slow passage through the resonance (with 
respect to both the electron transfer rate and the recip- 
rocal Rabi frequency) 



7rrj 2 e(r 



ML, 2 



MLA 



2\ji\y/& 



V 2 

1 M 



(62) 



Eq.© gives Q = e(T M L,2 ~ ^mla)^ > e at least for 
strong interaction when fl 2 3> r 2 ^ . 



Slow passage through the resonance and strongly chirped 
pulses 



Q = e- 



ML, 2 



MLA 



x exp 



-ttQ 2 

m 



Ias = 27re- 



ML.2 



r 



ML, l 



fi 2 



r M 

2 



(60) 



where W ia _i/ 2 (z) is the Whittaker function [39j. The 
graph of Ias as a function of Landau-Zener parameter 
and quenching parameter, which correspond to Q 2 / \ 
and I M j |/2 1, respectively, in terms of our representation, 
can be found in Fig.l of Ref . [38j . 

When chirp is fast with respect to the rate of the elec- 

tron transfer, -M- <C 1, one gets from Eq. 
I Ml 



ML, 2 



MLA 



r 



M 



1 — exp 



-nQ 2 



(61) 



where we have used the integral representation [39| of the 

1 1 2 

Whittaker function to calculate |lim 2 ^o W / ia,-i/2(^)| = 



Eq. (|62| can be obtained directly by integrating Eq. (l50l 
with respect to time for Lo(t) — loq — p,t and £1 — const 
Indeed, integrating Eq. (j50|) yields 



Q = 



I(t)dt 



M2 



Ml 



^M1^M2 
7rfi 2 e(r ML,2^ 'mR.,1 — ^ ML,!^ MR,2) 



IMI 



2 (rM2+r_Mi) 2 

Ml^M2 



+ (T MN /2f 



(63) 



In the special case T mi = ^ M2 = ^ m and B^{e 2 — £1) 
0, Eq.® leads to Eq.®. As a matter of fact, Eq.~ 
extends the case of slow passage through the resonance 
beyond the treatment of Ref. |38[ . 

Eq. (f6"3"| can be used for the excitation of a bridging 
molecule by Gaussian pulses, Eqs. (|57l) and (f58|) . as well 
when the pulses are strongly chirped 40j 



2|$»|»t 2 



(64) 
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For a strongly chirped pulse, one can ascribe to differ- 
ent instants of time the corresponding frequencies [ill ], 
i.e. different frequency components of the field are deter- 
mined via values of the instantaneous pulse frequency 
ui(t) for different instants of time. Then one can in- 
tegrate I(t)dt similar to Eq. ((63|) . bearing in mind 

that n(t) = (^)exp(-^ 2 t 2 ) « (^)eM-J^)^) 
is a much slower function of time than ui(t) = u>q — p,t rs 
u; Q + w ^t. Using Eqs.(j57l), |58|) . (f63]) and ([64]), we then 
get 

... Tm 2 + Tmi 

q « x 

Vtt(j-) 2 ^O c Ep T pO e (T M L,2^ M R,l — ^ML,l^MR,2) 
X — 

V^J ^"Ml 4r M ir M2 +\ 1 mn/1) 

(65) 

where £q « » since the magnitude 

X^o £ 2 {t)dt = 2/iQcEp — const is conserved in the course 
of chirping. Here is the pulse energy per unit area, /xo 
is the permeability constant, c - light velocity in vacuum. 
According to Eq. ([6"S"|) , in the case of slow passage through 
the resonance, Q ~ ( w )l f° r strong interaction 

when (#) 2 7if^1, (I >" + p Ml)2 » (r MJ v/2) 2 , and 
Q tends to a constant value for large |$"(o;)|. This 
elucidates the behavior observed in our simulations 
shown in Figs. 2] and [5] for moderate and large values of 
|$" (v)\. In addition, Eq. (|65]) explains why the growth of 
Q slows down for larger value of |$" (v)\ if pulse energy 
increases. 



VI. CONCLUSION 

In this work a theory for light-induced current by 
strong optical pulses in molecular-tunneling junctions 
have been developed. We have considered a molecu- 
lar bridge represented by its highest occupied and low- 
est unoccupied levels, HOMO and LUMO, respectively, 
and have derived a closed set of equations for electron 
populations of molecular states, molecular polarization 
and molecular excitation (exciton population) when two 
types of couplings between the molecule and the metal 
leads are presented: electron transfer that gives rise to 
net current in the biased junction and energy transfer be- 
tween the molecule and electron-hole excitations in the 
leads. 

We have used this formalism to analyze a novel con- 
trol mechanism by which the charge flow is enhanced by 
chirped pulses. For linear chirp and when the energy 
transfer between the molecule and electron-hole excita- 
tions in the leads is absent, this control model can be 
reduced to the Landau-Zener transition to a decaying 
level, which has an exact analytical solution. 

The relaxation parameters in the derived closed set of 
the equations of motion do not depend on the exciting 



electromagnetic field. This is true if the Rabi frequency 
ft is much smaller than the bath correlation frequency, 
uj c . If molecular states e m are far from the Fermi levels of 
both leads, uj c is determined by the frequency interval for 
the system-bath interaction matrix elements V^ K ^ and 

Vkk ,K ^ an d the density of states of metal leads. The last 
can be evaluated as 1-10 eV [2pj . As a matter of fact, the 
approximation of constant relaxation parameters, which 
do not depend on exciting electromagnetic radiation, is 
consistent to the RWA used in our theory. 

The situation is different if we assume that the molec- 
ular level position is pinned to the Fermi energy of a 
lead, that may lead to highly nonlinear current voltage 
dependence [16I | . In this case uj c is determined also by the 
frequency interval at which /if(e) is essentially changed 
that is ~ k B T/h (see Eqs.([2"D]) and (|2"g]l). In the last case 

can be of the same order of magnitude as oj c in the 
RWA, and the dependence of the relaxation parameters 
on exciting electromagnetic field [2(| must be included 
into the theory. 

To end this discussion we note that in this work we 
have investigated a model process driven by light absorp- 
tion in a molecular bridge connecting metal leads. As 
already discussed, the geometry considered is potentially 
advantageous because of the possible local field enhance- 
ment due to plasmon excitation in the leads. It should be 
emphasized however that other processes, not considered 
in this work, may play important roles in nanojunction 
response to incident light. First, direct electron-hole ex- 
citations of the metal leads [U, |42j may affect response 
in an adsorbed molecule that goes beyond the local field 
enhancement associated with plasmon excitation. Sec- 
ondly, experimental realization of strong local excitations 
in nanojunctions requires careful consideration of heating 
and heat dissipation and conduction (43[. Heating may 
be kept under control by driving the junction using a se- 
quence of well separated optical pulses, as envisioned in 
the proposed experiment, but it should be kept in mind 
that a more detailed consideration of this issue may be 
needed. 
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APPENDIX A 

Calculate the steady-state current in the absence of 
the radiative and nonradiative energy transfer couplings, 
Vp and Vn ■ The corresponding solution of Eq. ([55)) is as 
follows: n m — Wumj^ Mm- Substituting it into Eq. ([4"01) 
and using Eq. ([2~Tj) . we get 

I = e J2 TML f MR ^ [fR{e m ) - f L (s m )] (Al) 

~T» 1 Mm 

m— 1,2 

The last formula can be written 
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J = e > ^ !— / de[/i ? (e)-/ L (e)]J(e-e m ) = 

_e_ v- _ _ _ t^__ /" ^ f R {e)-f L {e) 



m=l,2 



(A2) 



using the well known representation for S(x) 



APPENDIX B 



S(x) = — lim 



7T x z + eH 



(A3) 



The limit limr Mm /2— >o on the right-hand side of Eq. (IA2p 
is consistent with the Markovian approximation in the 
sense that relaxation parameters rjf m /2 are small in 
comparison to the bath correlation frequency, u) c . 
The term limp Mm 

/2^0 ( £ -e m )!i fi -2 + [r Mr 011 the 

right-hand side of Eq. (|A2|) can be written 



77W 



lim 



1 



/2^o (e - e m ) 2 h- 2 + [T Mm /2] 2 



^"mm (^)^mm(^) 

(A4) 

where G^ m (£) and G^ lm (e) are the retarded and ad- 
vanced Green's functions, respectively, [l6| in the Marko- 
vian approximation. Substitution of Eq. (|A4[) into 
Eq. (|A2[) leads to the well known Landauer formula for 
the current 12111 . 



Let us show that n\ = 1 — n?, when Tmi — Tm2 = ^m- 
Summing Eqs. (|4Tj) and (|4^|) at given condition, we have 



dt 



(Bl) 



where we denoted y — n\ 
last equation is as follows 



n 2 



1. The solution of the 



y(t) = i/(0) exp(-r M t) 



(B2) 



where y(Q) = 0. This gives n\{t) + ri2(t) = 1 even in the 
presence of energy transfer when Bn(s2 — £i) ^= 0. 
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